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Abstract 

Using spin-density- functional theory, we study the electronic states of a two-dimensional 
parabolic quantum dot with up to iV = 58 electrons. We observe a shell structure for 
the filling of the dot with electrons. Hund's rule determines the spin configuration of 
the ground state, but only up to 22 electrons. At specific N, the ground state is de- 
generate, and a small elliptical deformation of the external potential induces a rotational 
charge- density- wave (CDW) state. Previously identified spin-density-wave (SDW) states 
are shown to be artifacts of broken spin symmetry in density-functional theory. 



PACS numbers: 73.23.-b, 73.23.Hk, 73.61.-r 



1 



Quantum dots have recently attracted much interest both experimentally and theoret- 
ically. One realization of a quantum dot is a small island fabricated in a two-dimensional 
electron gas laterally confined by an external potential and containing a few to a few hun- 
dred electrons 0. Experimentally, measuring the tunnel conductance^ and capacitance 
by changing the gate voltage attached to the quantum dot, one observes a peak every 
time the average number of electrons increases by one. The spacing of peaks, or addition 
spectrum, reflects the energy differences between ground states of the dot with different 
numbers of electrons. Each disordered dot has its own characteristic addition spectrum, 
but recently it has become possible to fabricate dots so clean that the addition spectra 
are reproducible from dot to dot||. Among the features of these clean, parabolic dots are 
atomic-like shell structures, Hund's rules, and reproducible transition rates @]. 

The advent of atomic-like spectra in quantum dots calls for appropriately quantitative 
theoretical tools. Presently, exact diagonalization of the full Hamiltonian is limited to a 
small number of electrons in the dot|5|, Thomas- Fermi, Hartree@, and Hartree-Fock 
methods [0, H all suffer from sizeable systematic errors. Here, we treat the electronic 



states of a dot using the spin- density-functional methodjlOl, [IT], [L2|]. We find shell struc- 
tures in the addition-energy spectrum for a circular, parabolic external potential. Hund's 
rule determines the ground-state spin configurations, but only up to 22 electrons. Ellip- 
tically deforming the external potential eliminates the shell structures, and Hund's rule 
is replaced by a more Pauli-like behavior of the total spin. At specific N, the ground 
state is degenerate, and a small elliptical deformation of the external potential induces a 
rotational charge- density- wave (CDW) state. The spin-density-wave (SDW) states found 
by Koskinen, Manninen, and Riemann (KMR) Jl2| are artifacts of broken spin symmetry 
in density-functional theory. 

We solve the following Kohn-Sham equations numerically for a two-dimensional parabolic 
quantum dot, and iterate until self-consistent solutions are obtained [p~3|l ; 
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Here a denotes the spin index, ((r) is the local spin polarization, and E xc is the exchange- 
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correlation energy functional, for which we use the local- density approximation fl4 



E xc = J p(r)e xc [p(r),C(r)]rfr (3) 

c(r) = w^gu. (4) 

To solve the equation, we expand the (r) in a Fock-Darwin representation; 
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guerre polynomial, I = \jTi/ {2m*ujQ) and Xc is a spin function. The non-interacting, 
single-particle levels form a ladder, = (2n + \k\ + l)ftu;o = Mfiu, with rung de- 
generacy M. The ground-state energy of a quantum dot with N electrons is obtained 
from 

We use the material constants for GaAs, m* = 0.067m, k = 12.9, and the external po- 
tential is fixed at Hujq = 3.0 meV. The resulting dimensionless interaction strength is 
(e 2 / kIq)/?iujq = 1.9, where £ = J%~/ (m*uj ). 

Shell structure - At low temperatures, electron hopping into a dot containing iV elec- 
trons is suppressed except when the ground-state energy E(N) is equal to E(N+l). This 
degeneracy condition determines the observed conductance oscillation peaks which occur 
at the chemical potentials fi(N + 1) = E(N + 1) — E(N). The addition energy A(JV) 
needed to put an extra electron in the dot is obtained from A (AT) = /i(iV + 1) — fi(N) = 
E(N + 1) - 2E(N) + E(N - 1). Fig. 1(a) shows the addition energy A(N) as a function 
of electron number N for a circular, parabolic potential. The dotted line indicates A (AT) 



obtained from a classical electrostatic analysis with no kinetic energy [15]. Overall, A(AT) 
decreases with N, as the dot and its capacitance grow. On average, the addition energy 
obtained from the density-functional calculation is close to the classical electrostatic re- 
sult e 2 /C. However, we see small zig-zag structures, and large peaks at electron numbers 
N = 2,6, 12, 20, 30, 42, and 56. In the single-particle spectrum for a parabolic potential, 
the electronic states of the dot form closed shell structures at these numbers. Even in the 
presence of electron-electron interaction, extra energy is required to add one more electron 
to a closed shell. The peak heights decrease as the number A^ increases, consistent with 
recent experiments which observed peaks in the addition energy up to 12 electrons [Q. 
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Hund's rule - By analogy to atoms, we expect that Hund's rule for total spin will apply 
in the present situation. According to Hund's rule, as degenerate states are filled, the total 
spin 5 takes the maximum value allowed by the exclusion principle and becomes zero for 
closed shells. Fig. 2(a) shows the spin configuration as a function of the electron number 
TV for the circular, parabolic potential. The dotted line represents the spin configuration 
when Hund's rule is satisfied. We can see that, for up to 22 electrons filling the dot, the 
spin configurations obey Hund's ruleQ. For larger dots, Hund's rule is violated and the 
high spin states are suppressed. In particularly, the total spin becomes zero at electron 
numbers N = 24,34,46 instead of the expected 5 = 2. In the upper inset of Fig.l, we 
show the total energy E(N) as a function of the total spin for these states. At N = 16, the 
5 = 2 state is 0.10 meV lower in energy than the 5 = state, which follows Hund's rule. 
In contrast, at iV = 24 the 5 = state is 0.05 meV lower in energy than the 5 = 2 state. 
This trend is enhanced as the number increases. These energy differences are sufficiently 
small that weak magnetic fields, of order 300 gauss, will favor an 5 = 2 ground state. 

The breakdown of Hund's rule is due to the non-parabolic effective potential caused 
by Coulomb interactions. For example, without interactions N = 24 corresponds to 20 
electrons in filled inner shells and four "valence" electrons distributed among 10 degenerate 
states: \n, k) = |0, ±4), |1, ±2), and |2, 0), spin up and down. Coulomb interactions deform 
the radial potential and lower the energy of the single-particle states with larger angular 
momentum \k\. The system could minimize its exchange energy by creating an 5 = 2 
state, i.e., putting all four valence electrons into spin up states. Instead, for iV = 24, 
the system minimizes its single-particle energy by putting all four electrons into k = ±4 
states, giving a total spin 5 = 0, and breaking Hund's rule. 

To confirm this result, we performed an exact diagonalization of iV = 4 valence elec- 
trons in a restricted basis set of eight states: |0, ±4) and |1, ±2), spin up and down. The 
Hamiltonian we employed is 



N=4 
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i=l 

where the 'jrf term is introduced to split the degenerate single-particle energies eo,±4 and 
ei != |-2- The resulting eigenstates of the four electrons can by labelled by total spin 5 and 
S z and total angular momentum L z . In the lower inset of Fig. 1, we have plotted the 
total energy as a function of A = e% ±2 ~ eo,±4- If the splitting A is small, the ground 
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state is S = 2,L Z = 0, consistent with Hund's rule. But for A larger than 1.4 meV, the 
ground state becomes S = 0, L z = 0, indicating a violation of Hund's rule. 

Spin-density-wave states - For a dot with circular symmetry, the eigenstates can always 
be chosen to have definite angular momentum L z , and hence circularly symmetric charge 
density. Nevertheless, KMR|12| reported recently on a spontaneous breaking of circular 
symmetry in a spin-density-functional calculation of a parabolic quantum dot. Indeed, 
we confirm that Eqs. (1-5) yield spin- density- wave (SDW) ground states at particular 
numbers of electrons, e.g. N = 24,34, as reported in Ref. ||12|| . These are precisely the 
5 = ground states discussed above in the context of breaking of Hund's rule. Within 
spin-density-functional theory, even for S = 0, the system lowers its exchange energy 
slightly by mixing in k = ±2 states with the lower lying k = ±4 orbitals. The result is an 
SDW state. However, from our exact diagonalization studies with N = 4 in the restricted 
basis set, we find that the SDW states are due to an unphysical mixture between states of 
different total spin: S — 0, S z = and S = 1, S z = 0. Hence, the SDW states are artifacts 
of the well known difficulty of spin- density-functional theory that only the S z component 
of total spin can be specified. We conclude that the correct ground states for N = 24, 34, 
and 46 have S = 0, L z = and retain circularly symmetry. 

Charge- density-wave states - We also find that for certain N (cf. Table 1), Eqs. (1-5) 
predict a rotational charge- density wave (CDW) near the edge of the dot. Fig. 3 shows 
an example of such a CDW state for N = 31. The numbers shown in bold face in Table 
1 indicate a strong spin-density modulation, as for p^(r) at iV = 31, while a weaker 
modulation occurs in the opposite spin density. The numbers in bold face correspond to 
a closed shell plus one electron, indicating that the extra electron is added to the lowest 
orbital in the next shell, namely the one with highest angular momentum. By circular 
symmetry, this orbital is doubly degenerate. Hence the ground state of the entire dot is 
doubly degenerate. (In contrast with atoms, the spin-orbit interaction in GaAs dots is too 
small to split this degeneracy [[TBI . ) The spin-density-functional result is a mixture of these 
two degenerate ground states. For example, at N=31 the total angular momentum will 
be L z = ±5, giving a charge density modulation ~ | exp(i56*) + exp(— i59)\ 2 ~ cos 2 (56') 
as observed in Fig. 3. 

We have investigated the charge density p a {v) for N = 3 by exact diagonalization 
to confirm the above interpretation. As expected, we find that there are two degenerate 
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ground states, with L z = ±1, and that a coherent mixture of these states produces almost 
exactly the same charge density obtained in the density-functional calculation. 

Elliptical dots - To investigate the effect of removing circular symmetry, we consider 
elliptically deformed potentials V(r) = m*{J 1 x x l + ujyy 2 )/2, with cu 2 = (u; 2 + which 
lift the large degeneracies of the shell structure. Fig. 1(b) and (c) show that as the defor- 
mation grows the regular zig-zag pattern found for the circular potential becomes irregular 
and the large peaks at large N disappear. One can still see large peaks at iV = 2, 6, 12, 
and 20 electrons in (b), which are the remnants of the closed-shell structures. However, 
in (c) such large peaks are present only at N = 2 and 6. 

Fig. 2(b) and (c) show the spin configurations for the same deformed external poten- 
tials. We can see that Hund's rule is satisfied up to N = 15 in (b) but only up to only 
N = 8 electrons in (c). The high spin states are suppressed as the deformation becomes 
significant - the loss of the closed-shell structures results in a more Pauli-like behavior of 
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Fig. 3 shows the up-spin densities in deformed potentials for N = 31, = 16. We 
find that true CDW ground states are induced by the increasing elliptical deformation 
of the external potential. The ground-state spin is S = 1/2 in all three cases. The 
charge- density wave has period ~ cos 2 (58) and results from the mixing of the degenerate 
L z = ±5 states by the elliptical external potential. 

In conclusion, we have studied the electronic states of quantum dots with up to 58 elec- 
trons for parabolic circular and elliptical external potentials, using spin-density-functional 
theory and exact diagonalization. For a circular potential, we observe a shell structure for 
the filling of the dot with electrons. Hund's rule determines the spin configuration of the 
ground state up to 22 electrons. For specific numbers of electrons, CDW states appear 
on small elliptical deformation of the external potential, while previously identified SDW 
states |12[ are found to be artifacts of broken spin symmetry in density-functional theory. 
For elliptical potentials, the shell structures are lost with increasing deformation, and the 
spin configurations change from Hund's rule to a more Pauli-like behavior. 

We acknowledge O. Again, I.L.Aleiner, B.L.Altshuler, D.J.Chadi, Y.Meir, W.Kohn, 
and M.Stopa for comments and suggestions. 
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Figure Captions 

Figure 1: Addition energy A (AT) as a function of electron number N in the dot for (a) 
a confining parabolic potential nfuj^r 2 /2 and (b,c) elliptical confining potentials 
V(r) = m*{uj 2 x 2 + u 2 y 2 )/2. The dotted lines indicate the addition energy accord- 
ing to a classical electrostatic analysis |L5j. The parameters are (a) huo = 3.0 meV, 
(b) Uy/u) 2 = 11/13 and (c) uJy/uJ 2 = 5/7. (b) and (c) are shifted by 1.0 meV and 
2.0 meV, respectively. 

Upper inset - total energy in meV as a function of total spin for electron numbers 
N = 16, 24, 34, and 46 in (a). The origin of energy for each N is arbitrary. 
Lower inset - total energy for N = 4 electrons obtained by exact diagonalization 
within a restricted Hilbert space, |0, ±4) and |1, ±2), as a function of single-particle 
level splitting A. The energies E(S = 0, L z = 0) and E(S = 1, L z = ±2) are plotted 
relative to E(S = 2, L z = 0). 

Figure 2: Ground-state spin as a function of electron number N for (a) the parabolic confining 
potential and (b,c) the elliptical confining potentials of Fig. 1 The dotted line in 
(a) indicates the spin configuration when Hund's rule is satisfied. 

Figure 3: Charge-density distributions at iV = 31. Left column - from top to bottom, spin- 
up (iV^ = 16), spin-down (N^- = 15), and total charge-density distribution. Right 
column - spin-up charge-density distributions (N^ = 16) for elliptical potentials, 
UJ 1/ ljJ 'x = H/13 (top), ujy/u 2 = 4/5 (middle), and to 2 /col = 5/7 (bottom), respec- 
tively. 
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Table 1 
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